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Abstract 



Inflationary cosmologies predict Gaussian primordial fluctuations, but 
subsequent gravitational lensing of the CMB disturbs its Gaussianity. Knowl- 
edge of the specific signature of lensing is necessary to distinguish a lensed 
Gaussian sky from an intrinsically non-Gaussian one. In this paper we inves- 
tigate the pattern of non-Gaussian features in the CMB resulting from grav- 
itational lensing, both analytically and numerically. We describe the lensed 
CMB temperature distribution using the formalism of generating functionals 
and obtain the complete picture of emerging non-Gaussian effects. To test 
for the signature of lensing in practice, we propose a scale-dependent statistic 
based on cumulants of the CMB temperature, and we compute that statistic 
numerically for typical inflationary models. The effect is expected to be most 
significant at angular scales ~ 10' and its magnitude is within observational 

capabilities of future satellite missions. 
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I. INTRODUCTION 



The cosmic microwave background (CMB) radiation is a direct imprint of the primordial 
matter, and its importance for cosmology cannot be overestimated [|IJ. The forthcoming 
MAP 0j and PLANCK || satellite missions are expected to carry out high-resolution, 
low-noise measurements of the CMB temperature and polarization, yielding a stringent 
test of our theories of cosmological structure formation. One of the pronounced differences 
between the competing theories is a Gaussian distribution of the density fluctuations in most 
inflationary theories versus the generically non-Gaussian inhomogeneities in scenarios based 
on topological defects. For this reason, detection of non-Gaussian features in the CMB is 
important for distinguishing the cosmological models. 

However, even if the fluctuations of the CMB were Gaussian at the surface of last scat- 
tering, their later evolution would have introduced certain non-Gaussian features in their 
distribution. One of such effects is gravitational lensing by matter present along trajecto- 
ries of the primordial photons. While the effect of weak gravitational lensing on the power 
spectrum of the CMB fluctuations has been extensively studied (see e.g. |4| and references 
therein), including its relevance to observations of large-scale structure and cosmological 
parameter estimation || , it is only recently that attention has been drawn to the appearance 
of non-Gaussian signal in the CMB due to weak lensing J7]||. 

To identify a sky map as coming from an intrinsically non-Gaussian model rather than 
from a lensed inflationary one, we need to distinguish the non-Gaussian signature of lensing 
from that of other processes. The goal of this paper is to investigate in detail the effect of 
weak gravitational lensing on the distribution of the CMB fluctuations and to fully charac- 
terize the emergence of non-Gaussian features. Rather than restrict ourselves to a specific 
criterion of Gaussianity at the outset, we describe the lensed sky in a general way by using 
the formalism of generating functionals, somewhat similar to that used in quantum field 
theory. A description of a random field in terms of a generating functional is equivalent 
to giving the full probability distribution and in principle is enough to evaluate any field 
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statistic. 

To test for the signature of lensing in practice, we chose a criterion based on cumulants of 
the temperature distribution ||. Cumulants are generic additive and scale-sensitive probes 
of Gaussianity defined through moments of the field and are easy to compute in our con- 
text. Recent studies of geometrical criteria of Gaussianity [TIJ show that the "generalized 



genus" statistics based on Minkowski functionals [|llj] are sensitive to the 3-point temper- 
ature distribution at points separated by the smallest map scale. A direct examination of 
that temperature distribution using cumulants yields a statistic essentially equivalent to 
or possibly more discriminative than the generalized genus. Thus motivated, we apply a 
cumulant statistic to detect the non-Gaussian effects of lensing. 

The distribution of the CMB after lensing depends on the initial CMB power spectrum 
and on the distribution and evolution of the lensing sources, all of which have to be taken 
from a particular cosmological model. We compute the effects of lensing in typical infla- 
tionary CDM cosmological models, using the program CMBFAST |L2| to obtain the unlensed 
CMB power spectrum and the lens distribution. As our Gaussianity test, we evaluate the 
dimensionless cumulant X22{9) which of all cumulants is expected to show the largest devi- 
ation from the Gaussian value Xij = 0> as we U as the next largest cumulant x.33- Using a 
wide range of cosmological parameters, we find that the maximum signal as measured by 
X22 {&) is generically expected to be of order 10 _2 -10 -3 at angular scales 9 ~ 10-20 arcmin. 
The value of X33 is smaller by several orders of magnitude. 

The paper is laid out as follows: In the next section we present an outline and the 
results of our main calculation of the generating functional for the temperature distribution 
after lensing, in a model with Gaussian initial fluctuations and Gaussian distribution of 
lensing sources. The detailed definitions of the formalism are found in Appendix [A], while 
the calculation itself is deferred to Appendix JTJ. We also comment on a generalization of 
our formalism to non-Gaussian initial conditions and lens distributions. In Section [TTT] and 
Appendix |Cj we compute the observable effect of lensing as measured by the cumulant X22 (#) 
and investigate the possibility to detect this signal, which may just be at the lower limit of 



angular resolution and sensitivity of the forthcoming PLANCK mission. 



II. NON-GAUSSIAN EFFECTS OF LENSING ON THE CMB 

To achieve a clear understanding of non-Gaussianity due to lensing, we first wish to 
concentrate on the simplest case when both the initial CMB fluctuations and the lensing 
sources are Gaussian. The effect of lensing is to make a primordial photon appear to originate 
from a slightly different direction of the sky. Accordingly, one uses the "lensing field" g (x) 
to shift the original temperature field T (x), so that after lensing it becomes the observed 

r(x), 

T(x)=T (x + g(x)). (1) 

The temperature before lensing is taken as a homogeneous Gaussian random field T (x) 
with zero mean and known correlation function 

(T (x)T (y)> = Co(x-y). (2) 

(Here x is a point in a region of the 2-dimensional Euclidean space R 2 which represents 
a small flat patch of the sky. As we shall see below, all lensing effects appear on sub- 
degree scales, and the small-angle approximation used throughout this paper suffices for 
our purposes.) The "lensing field" g (x) is also assumed to be a homogeneous zero- mean 
Gaussian random function on R 2 with values in the same space R 2 and known (matrix- 
valued) correlation function 

C g (x - y) = C« (x - y) = i*«C, (x - y) . (3) 

(The function C g is the same as C g i in the notation of Ref. Q.) Due to homogeneity, all our 
correlation functions depend only on distance, C (x) = C (|x|). 

We shall assume below that the angular displacement g due to lensing is "small" and 
it makes sense to expand in g. Quantitatively, the lensing effect is small if the standard 
deviation of the deflection angle g gives a scale at which the original temperature field does 
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not change appreciably. This condition can be expressed through the real-space correlation 
functions of the two fields as 

c s( o,«^|L. (4, 

As we verified numerically, this condition holds for all inflationary models. (The values of 
C g (0) are < 10~ 7 while the right hand side of Eq. (f|) is typically ~ 1CT 5 .) 

The correlation functions Cq (x) and C g (x) can be computed from any given cosmological 
model. The CMB correlation function is conventionally expressed through the angular power 
spectrum Cj, 

Co {0) = —Y, (2/ + 1) QPi (cos6) . (5) 
z=o 

The angular displacement g (x) is found by integrating the transverse gradient of the grav- 
itational potential (f) from the time of last scattering to the present time. The gravitational 
potential is related to the fluctuation of the primordial matter density 5p by a Poisson 
equation. The explicit expressions are derived by several authors (see e.g. |!3iH) an d we 
shall not need to reproduce them here. Although the distribution of density fluctuations Sp 
and hence of the gravitational potential 4> may be non-Gaussian due to either primordial 
topological defects or nonlinear evolution, the averaging over photon trajectories is expected 
to wash out, due to the central limit theorem, at least some of the non-Gaussianity. We 
shall comment on more general initial conditions later, and for now we shall assume that a 
certain model of cosmological structure formation has been chosen and the corresponding 
correlation matrix C 9 of g (x) has been computed. 

Even if both the original field To (x) and the deflection g are Gaussian, the field T (x) 
after lensing has a non-Gaussian distribution. We shall describe the distribution of T (x) 
through its generating functional Z [J] defined as the average 



Z[J\ = ( exp 



J fx) T fx) rfx 



= (exp [-U X T X ]) T , (6) 

T(x) 



which depends on the parameter function J(x). The generating functional of a Gaussian 
distribution must be a Gaussian functional, such as 



Z G [J] = exp 



J(x)J(y)C(x,y)dxdy 



exp 



T T (~ixy 

-0 X Oy\-J 



(7) 



A distribution with slight deviations from Gaussianity may be described by a generating 
functional of the form 

/ -\3 



Z [J] = exp 



1 

2 "x" 



z. t t c x y 



3! 



- J X JyJ + 



(8) 



where the Gaussian (quadratic in J) term in the exponential still gives the correlation 
function C (x — y), while the coefficients C n , n = 3, 4, ... of the non-Gaussian terms are the 
general cumulants (see Appendix [A| for a more detailed exposition). Taken together, the 
cumulants fully characterize the deviations from a Gaussian distribution, and to find them 
or equivalently the functional of Eq. (|8D is the main goal of this section. 

We now give a sketch of the calculation of Z [J] for the lensed temperature field; the full 
details are found in Appendix |B|. Combining Eqs. (0) and (|6]) gives 



Z[J] 



exp 



J(x)T n (x 



(x)) <ix 



(9) 



here the average must be performed over all realizations of 7o(x) and g(x). The result can be 
expressed through the unmodified correlation function Co (x) and the generating functional 

1. 



exp 



2 -H^CfH, 



(10) 



of the deflection field g, and we obtain (cf. Eq. ( |E7| 



'^H(y) ) Jy 



Zg [H] 



(11) 



H=0 



We are actually more interested in finding the logarithm of Z [J] since, as Eq. (|§D illustrates, 
In Z [J] is the generating function of the cumulants. If one expands Co in the exponential of 
Eq. (P), 



Z [J] = exp 



1 

2' 



t r xy t — 



iJ x J y VC (yL-y) 
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J x J y VVC (x - y) 
2 • 2! 



5 



5H(x) <JH(y) 

2 



5H(x) 5H(y) 
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+ ... 



Z 9 [H] , 



(12) 



H=0 



one obtains an expression formally similar to the generating functional of Green's functions 
in interacting quantum field theory with the "couplings" represented by the coefficients at 
the functional derivative terms. Following the analogy with quantum field theory, we notice 
that In Z [J] is the sum of the appropriate connected Feynman diagrams having no external 
lines (as shown in Figs. 1,2), with vertices being derivatives of Co (x) of all orders and the 
"propagator" being C g (x — y). The full expression corresponding to a vertex of arbitrary 
order is given by Eq. (|B9|) . Unlike the usual perturbation theory, the order of a diagram is 
now equal to the number of propagators in it, because the smallness of the lensing effect is 
passed on to the propagator rather than to the vertices. Since the expansion in Eq. fljjD is 
formally to all orders, an infinite number of vertices with arbitrary number of legs arises, but 
diagrams with many-legged vertices are suppressed because they require more propagators 
and therefore are of higher order in C g . 

Since all vertices in our diagrams are quadratic in J, each diagram yields an even power 
of J. After adding up all diagrams, we therefore expect to obtain an expression for the 
generating functional such as 



i 2 , , „„,„. i A , , , , 



In Z [J] — 2j JxJyC V + ^ J wJ xJ yJ zC ^ ^ 4~ J u J v J w J x JyJ z CQ ... (13) 

containing only even powers of J, with the even-order cumulants C, C 4 , C 6 , ... found as 
sums of diagrams with the appropriate power of J. The new correlation function C (x) and 
the fourth-order cumulant C4 are (to lowest order in C g ) given by the two diagrams in Fig. 
lb, 

C (x) = Co (s) + \ (C g (0) - C g (*)) (c' ' (x) + + O {C 2 g ) , (14) 

C A (w, x, y, z) = -C' Q (|w - x|) C (|y - z|) C g (|w - y|) (w " X) ' (y " Z) 



4 UVI " UVI " 17 yvi uu |w-x||y-z| 

+ (23 permutations of w, x, y, z) + O (C 2 ) . (15) 

(The expression in Eq. flT5p is effectively symmetrized in all four parameters.) The higher- 
order cumulants are more complicated and we shall not write them out explicitly. The 



functional of Eq. (|i~3|) together with the diagrammatic method of computing the coefficients 
C, C 4 , C 6 , ... is our main result describing the characteristic non-Gaussian contributions of 
gravitational lensing. 

From the spatial dependence of the terms in Eq. ([15]) it is clear that the 4-th order cumu- 
lant C4 vanishes for certain combinations of points, for instance C4 (x, x, x, y) = 0. It can 
be shown that similar restrictions hold for higher-order cumulants and, in particular, that 
the one-point cumulants C n (x, x, ...,x) for n > 2 vanish, which means that the one-point 
temperature distribution remains Gaussian. This is a manifestation of the specific signa- 
ture of gravitational lensing: it only modifies the distribution of relative angular distances 
between the CMB photons. 

The condition that first-order corrections to the power spectrum due to lensing are small 
is, by inspection of Eq. (|14"D, equivalent to Eq. (||), so the expansion in C g is justified. 
In typical cosmological models, the magnitude of C g is never larger than 10 -2 in units of 
the right hand side of Eq. (|4]). Since the total power of J in a diagram is equal to twice 
the number of vertices, it follows that diagrams contributing to higher-order cumulants 
will be suppressed by additional factors of C g , and the non- Gaussian signal in higher-order 
cumulants rapidly becomes negligible. The most significant effects of gravitational lensing 
are therefore expected in the power spectrum (i.e. in C (x)) and in the 4-th order cumulant 
C A . 

Eqs. (p^)-([T5|) are consistent with the results of Seljak [[| and Bernardeau ||. Our 
approach is somewhat similar to that of Ref. || where the non-Gaussian contributions to 
the 2- and 4-point correlators of temperature were investigated by expansion in the angular 
displacement g, and the results of Ref. || are of course reproduced by the leading-order 
expressions given by our Eqs. (|1^)-(|T5|). As also noted by Bernardeau ||, the correction 
to the power spectrum given by Eq. ([14]) corresponds, in the lowest order in C g , to a more 
precise formula derived in Ref. Q. 

Having treated the simplest case of Gaussian fields, we now consider the possible gener- 
alization of our method to non-Gaussian initial conditions and lensing source distributions. 
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Such a generalization can be implemented by introducing new kinds of diagrams contributing 
to the cumulant expansion. The non-Gaussian effects may not be limited to even-numbered 
cumulants in that case. The generating functional of the lensed temperature can still be ob- 
tained from Eq. (|12|) , where one needs to put in a non-Gaussian generating functional Z g [H] 
to represent the lensing sources as well as a non-Gaussian term (e.g. cubic or quartic in J) 
in the exponential for each non- vanishing cumulant of the initial CMB distribution. Then a 
reasoning similar to the one presented here and in Appendix |B| will again lead to a diagram- 
matic representation of cumulants. If the initial CMB distribution has nonzero cumulants of 
odd order, vertices giving odd powers of J will also be present in the diagrams. On the other 
hand, non-Gaussianities in the lensing sources lead to many-legged "propagators" (Fig. 3). 
For a generic weakly non-Gaussian distribution, one would expect the magnitude of an n-th 
order cumulant Cg of the lensing sources to be of order (C g ) n ^ 2 . Then, despite the infinite 
proliferation of vertices and propagators, the expression for any given cumulant will always 
be finite in all orders in C g and the leading order term will usually be given by few diagrams 
containing low-order vertices and propagators. 

Although the signature of lensing becomes less clear-cut for a non-Gaussian initial CMB 
distribution, one can still distinguish some intrinsic non-Gaussian features from those due 
to lensing. For instance, one can show that the cumulant C n (x, x, ...,x) is unchanged after 
lensing. 

Our formalism seems to require knowledge of the full generating functionals Z [J] and 
Z g [H] for the lensing sources and original CMB distributions, but in fact leading-order 
results can be obtained using only the first few cumulants of those distributions, such as the 

(3) 

three-point cumulant of the lensing sources C g which contributes to C4 according to the 
diagram in Fig. 3b. The cumulants of the lensing sources can in principle be estimated from 
iV-body cosmological simulations or from galaxy surveys, which is an interesting direction 
for future work. 
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III. EXPECTED EFFECT IN GAUSSIAN CDM MODELS 



Many different statistics have been proposed and used to probe for non-Gaussian features 
in the CMB. Here we shall consider a form of the cumulants-based criterions suggested in 
which in a way generalizes the skewness and kurtosis [H] and the genus statistics |]T5| . 



Our particular form of the cumulant criterion is related to the morphological descriptors 
of Gaussianity based on Minkowski functionals 0,0, which are a generalization of the 
topological genus. Minkowski functionals have been applied in Ref. [TO to regularly pixelized 
maps, and it was found that all three morphological descriptors (the area, the boundary 
length, and the genus of temperature excursion sets) are effectively determined by the three- 
point distribution of temperature at nearest neighbor pixels in the map. One might as 
well inspect this three-point distribution directly and investigate whether it is Gaussian; 
such analysis would be in principle equivalent to inspection of the Minkowski functional 
curves calculated from the same map. A standard method of checking that a distribution is 



Gaussian is by computing its cumulants [16j . Thus we arrive to a criterion of Gaussianity of a 
CMB map which consists of examination of various cumulants of the three-point distribution 
of temperature at neighbor pixels. Note that the evaluation of a given cumulant requires to 
compute certain moments of temperature at neighbor pixel points, and since there is a fixed 
number of neighbors for each pixel, the total number of necessary operations is linear in the 
number of pixels in the map. 

For the purposes of this paper, however, we shall limit ourselves to considering the two- 
point cumulants of the distribution of temperature, computed for pairs of points separated 
by a fixed angular distance 9. The cumulants of a two- variable distribution are labeled by 
two indices and may be defined by 



~X,mn (^) i 



m+n 



lnf e (p,q) (16) 

x,y=0 



dp m dq n 

where fg (p, q) is the generating function for the joint distribution. To be definite, we denote 
the two field values by / (0) and / (6), although we imply in our expressions below that 
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averaging is performed over all pairs of points separated by 9. The formulae for the practical 
evaluation of cumulants from a given map follow from Eq. (0); they are simplified if the 
field / (x) = AT (x) /T has zero mean, as we shall assume. The cumulant xu corresponds 
to the lensed power spectrum, while Xmn with m + n > 2 represent the non-Gaussian signal. 
We shall be particularly interested in the cumulant 

X22 (9) = (f (0) 2 / (6) 2 ) - 2 (f (0) / (6)} 2 - (f (0) 2 > (/ (9) 2 ) (17) 

which characterizes the (appropriately collapsed) 4-point correlation function of tempera- 
ture. Since the cumulants defined by Eq. (|T6|) have dimension f m + n j it is convenient to 
normalize them to units of standard deviation <Tq = (/ (0) 2 ), 

X mn (0) = X mn(e)a^ m+n) . (18) 

The expected value of the cumulants Xmn (9) for the lensed sky is computed in Appendix 
|Cj where it is shown that the only nonzero cumulants are those with m = n. To characterize 
the magnitude of the effect it is convenient to introduce the parameter 

a 2 (9) = C g (0) - C g (6) (19) 

which describes the relative deviation of photons initially separated by angular distance 9 
(it corresponds to cr 2 (9) of Ref. Q). Up to O (o~g), the first few cumulants of the lensed sky 
are (Eqs. (|U8D-( PT0|) ) 

Xu (9) = C (9) + a 2 



+ 0(a e g ), (21) 

X33W = 72^(Q 2 ^' + 0(a 9 6 ). (22) 

(In these expressions we imply that a g and all derivatives of the correlation function Cq (x) 
are evaluated at x = 9.) Cumulants Xnn for n > 3 are of order O (o^) or smaller, while all 
other cumulants vanish. This is the specific signature of the gravitational lensing. 
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r" 4- -° 
Go+ e 



+ o( 



(20) 



X22 (9) = 4a 2 (Q 2 + 4a 



Now we turn to numerical estimates of the non-Gaussian signal measured by X22 (#)• Our 
objective is to find the angular scale 9 at which the effect is most significant and to determine 
whether the signal is within our observational capability. To obtain a numerical answer, we 
need to choose the parameters of a particular cosmological model (the most obvious ones 
to vary are h, Q^, Q\ and n s ) and compute the two correlators, Cq (9) for the CMB before 
lensing and a g (9) for the lensing sources. The resulting functions are then substituted 
into Eqs. (p0|)- (|22]) . The CMB power spectrum 1(1 + 1) Gij (2tc) can be generated by the 
program CMBFAST |l2j and converted to C (9) using Eq. (^); the same power spectrum is used 



to estimate the Gaussian cumulant variances (cf. Eq. (|C17|) ). The latest version of CMBFAST 
also computes the effect of gravitational lensing on the power spectrum following Ref. [[|; 
the author used the CMBFAST code to directly output the lens correlation function a g (9). 

The resulting dimensionless cumulants X22 (9) and X33 (9) are plotted in Fig. 4. We chose 
two cosmological models which show the typical range and behavior of the cumulants. The 
generic prediction for all models is a peak in X22 of magnitude ~ 10 _2 -10~ 3 at ~ 10'-20 / . 
In open models (Q < 1), the peak is generally lower and occurs at larger observed angular 
scales than in flat models. The numerical results also show that the higher cumulants (as 
well as the higher-order terms in X22) are several orders of magnitude smaller, as expected. 
The variance of the cumulant X22 is typically ~ 10~ 3 , which as our graphs show only allows 
detection at scales near the angle 9 — 9* of the maximum signal. It is necessary to note 
that the variance of the cumulant X22 given by Eq. (|C17| ) is computed for the theoretical 
limit of infinite-resolution, noiseless, full-sky maps, and a real experiment will give a larger 
error bar. If the temperature points were completely uncorrelated and pixels were squares 
of angular size 9q, the distribution of pairs at angle 9 would be estimated from N ~ 8tc9*/9q 
independent pairs. For the relative standard error 1/yN to be less than 10 -3 , the map 
resolution must be 9 < 15'. This angular scale is beyond the capabilities of MAP but is at 
the lower limit of the PLANCK resolution. 

It is also easy to see that pixel noise present in the data does not affect the results as long 
as the signal-to- noise ratio is somewhat larger than 1. Since the cumulants are additive, the 
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observed value of X22 is the sum of the cumulant values of signal and noise. Uncorrelated 
pixel noise will only contribute to the variance of X22 which depends on the fourth power of 
the correlation function. Therefore the contribution to the variance due to pixel noise will 
be smaller than the intrinsic variance by the signal-to-noise ratio to the fourth power. 

In conclusion, we note that the non-Gaussian effect of gravitational lensing is a necessary 
component of the CMB signal in all cosmological models, and in inflationary scenarios it 
leads to deviations from the Gaussian distribution which are small but observable by next- 
generation satellite missions. The formalism of generating functionals uncovers the signature 
of gravitational lensing, and we have shown that certain cumulant-based statistics can detect 
it on subdegree scales. 
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APPENDIX A: WEAKLY NON-GAUSSIAN RANDOM FIELDS 

Here we describe the formalism of generating functionals we use to characterize the 
non-Gaussian features of the CMB. We also give a definition of cumulants. 

Recall that the generating function f(ki,...,k n ) of the joint probability density 
p (xi, x n ) for a set of random variables x = {xi, x n } is defined by 

/ (k) = (e- kx ) x = J exp {-iV Xl ) p (x) dx. (Al) 

The expectation value of any function A of x (treated as a formal series expansion in x) can 
be then obtained by differentiation, 
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.d_~ 
'9k 



(A2) 



k=0 



A description of the distribution of variables Xi in terms of the generating function / (k) is 
equivalent to the description in terms of the joint probability density. 

In the infinite-dimensional case of random functions ("fields") on a Euclidean space, we 
consider the generating functional Z [J] defined as the average 



Z [ J] = ( exp 



—i J J(x)/(x)<ix 



(exp {-U x f x }) 



/(x) 



/ 



(A3) 



over all realizations of /(x). If the generating functional Z [J] for a field /(x) is given, one 
can obtain the correlation functions of f(x) and averages of any functional A, 



(/(x 1 ).../(x B )>=*' 



5J(xi) ...<JJ(xi) 



Z[J) 



j=o 



(A\f {*)])= A 



S 



8 J x 



Z[J]. 



(A4) 



(A5) 



J j=o 



This implies the normalization (1) = Z [0] = 1. 

Let us give an example of how the generating functional Z [J] of a random field /(x) can 
be directly used to evaluate some field statistics, such as the ones used for CMB analysis. 
If one is interested in finding the distribution of values /(x) at certain given points xi, 
x n , one can obtain such distributions directly from Z [J]. Again, we should look for the 
generating function / (ki, k n ) for the distribution of values /(xi), /(x n ), as defined by 
Eq. (|A1[) , rather than for the probability density itself. It is easy to show that the generating 
function / (k\, k n ) is found by substituting into Z [J] a suitable function J(x): 



/ (kt, k n ) = Z [J(x) = (x - xi) + ... + k n 5 (x - x n )] 



(A6) 



Another example: The generating function / (fc , k) for the joint distribution of / and its 
gradient V/ taken at a reference point x is 



/ (fc , k) = Z [k 8 (x - x ) - k • V8 (x - xo)] • 



(A7) 
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We turn now to the question of characterizing non-Gaussian features of random fields. It 
is straightforward to show that the generating functional of a Gaussian distributed random 
field is a Gaussian functional, 



Z G [J]= exp 



-i J n(x)J(x)dx.- ~ J J(x)J(y)C(x,y)dxrfy 



exp 



IjlJ —J X C y Jy 



(A8) 



where /x is the mean value and C is the correlation function. If the field in question is 
homogeneous and isotropic, its mean value is a constant /i (x) = fi and the correlation 
function depends only on the distance, C (x, y) = C (|x — y|). A general non-Gaussian field 
will not have the generating functional Z[J] of the form of Eq. ( |A8| ) ; however, we could 
attempt to express Z [J] as an exponential of a formal series in J, written in condensed 
notation as 

\3 



Z [J] = exp 



-i) 



-J x JyC ^ ~\- JxJyJzC^ -\- 



3! 



(A9) 



with some functions C n (x 1; x n ), n = 3, 4, ... as parameters. For this representation to be 
possible, it is necessary for all moments of the distribution to be finite, which of course is 
not the case in general. However, we are interested in distributions which are "weakly non- 
Gaussian", and in particular, we expect all moments of such distributions to exist. So we 
shall take the possibility of a series representation (Eq. (|AD|)) for the generating functional 
of a random field as a definition of its weak non-Gaussianity. 

We regard the coefficients C n (x 1; ...,x n ), n — 3, 4, ... (which are in fact the cumulants 
of the distribution of /) as general descriptors of the non-Gaussian features of the random 
field /. The cumulants can be related to the correlation functions by rewriting Eq. 



as 



C n (xi, x n ) 



\nZ[J] 



(A10) 



j=o 



<5J(x 1 )...<5J(x 1 ) 

and replacing the resulting derivatives of Z [ J] at J = by the corresponding correlation 
functions. It follows that C\ = and C2 (x, y) = C (x — y) are the mean and the correlation 
function of the distribution, and for example the third-order cumulant is 
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C 3 (x, y, z) = (/ (x) / (y) / (z)> - [C (x - y) + C (y - z) + C (z - x)] /x + 2/i 3 . (All) 

Because In Z [J] is the generating functional of cumulants, similarly to the generating func- 
tional of connected Green's functions in quantum field theory, the cumulants are sometimes 
called the "connected moments" . For a Gaussian distribution, all cumulants C n with n > 3 
vanish by construction. 

A convenient property of the cumulants C n (to which they owe their name) is additivity: 
if / is a sum of two statistically independent distributions g and h, then the generating 
functional Zf for / will be the product Zf [J] = Z g [J] Z h [J] of the generating functional 
for g and h. Because of the logarithm in Eq. ( |A10| ), the cumulants for g and h will add to 
give the cumulants for /. 

The description of a non-Gaussian field in terms of the cumulants C n is a compact way 
to put together the information from all n-point correlation functions. 



APPENDIX B: GENERATING FUNCTIONAL FOR COMPOSITION OF FIELDS 



Here we compute the generating functional for the lensed field of Eq. ([I]). The average 
over the space of realizations of T(x) means the average over To (x) and g(x). We assume 
that both Tq(x) and g(x) are Gaussian distributed with generating functionals 



Z Q [R] = exp 



Z g [H] = exp 



i?(x)Co(x-y)i?(y)dxdy 



-~ J H T (x)C 9 (x - y) H(y)rfxdy 



= exp 



= exp 



^R X C Q Ry 



(Bl) 



vB2) 



(in the last expression, H is an i? n -valued function and matrix multiplication is implied 
between C g and H). 

According to Eq. (|A~5|), the generating functional for T (x) is formally written by substi- 
tuting i5/5R instead of T and iS/5H instead of g in Eq. 



Z[J] ^exp f-i j J(x)T (x + g(x))rfx^ 



T ,S 
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To make the calculation clearer, we first evaluate the average only with respect to T , 
which corresponds to taking the derivative 5 /SR. We temporarily introduce the inverse shift 
function g (x) which satisfies g (x) = g [x + g (x)] so that we can change variables 

J J (x) T (x + g (x)) rfx = j J (z - g (z)) T (z) dz. (B4) 

(We shall not need an explicit expression for g (x) in what follows.) Then we obtain 



exp(-* / J(x)T (x + g(x))dx) ) =Z [J(x-g(x))] 

To 



cxp 



rfx / dyJ(x-g(x))C (x-y) J(y-g(y)) 



cxp 



rfx / rfyJ(x)C (x-y + g(x)-g(y)) J(y) 



exp 



--J x c*y[ g ]j y 



(B5) 



The correlation function can be expanded as 

&* [g] = Co [x - y + g (x) - g (y)] = C7 + [g k (x) - g k (y)] V fc C (x - y) 

\g k (x) - g k (y)l \g l (x) - g l (y)l 
+ l V } 1 ;J L 1 ^ 1 JJ V fe VA (x - y) + O (g 3 ) . 



(B6) 



We now only need to substitute i5/5H instead of g to express the generating functional as 



Z [J] = exp 



-g-JaA ( x - y + '- 



5 



5 





exp 






H=Q 





(B7) 



5H(x) 5H(y), 

The expression of Eq. (|B7|) is of course understood in the sense of perturbative expansion 
as shown by Eq. (|B6|) , yielding a sum of terms with functional derivatives of increasing 
order. In the context of weak lensing, the correlation function C g can be regarded as a small 
parameter in the sense of Eq. and expansion in powers of C g is a natural choice. Each 
pair of functional differentiations 5/5H applied to the Gaussian generating functional Z g [H] 
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would add at least one more power of C g , therefore the expansion in C g is consistent with 
truncating the operator series at a finite order of functional derivatives. 

Exploring the quantum field theory analogy further, we regard Eq. ( |B~7| ) as the generating 
functional of Green's functions for an interacting quantum field, with an "interaction vertex" 
for each term of the expansion in gradients of Cq. The connected Green's functions of 
this theory are generated by In Z [J] , which in our context is the generating functional of 
cumulants. According to standard rules, we can express In Z [J] as the sum of all connected 
diagrams made up of vertices with any number of legs n > 1; the propagator is the 
matrix C g (x — y). Since all functional derivatives are taken at H = 0, the resulting diagrams 
are all "vacuum diagrams" without external legs (Figs. 1, 2). 

To obtain the general expression V^ n ' for the vertex with n legs, we need to consider the 
n-th term in the expansion of the correlation function, 

C (x - y + g ( X ) - g (y)) = C (x - y) + (g k (x) - g k (y)) V fc C (x - y) + ... 

+ ^ / Vfc i- Vfc » cr ° ( x - y) n p - x ) - 6 ( x * - y)i ny* w ^ + - ^ 

i i 

The required vertex coefficient is the factor at JX^*' ( x «) without the n! denominator, 
integrated with dxdyJ(x) J(y)/2, 

V£}. kn (x lf x n ) = I J dxdyJ(x) J(y)V fcl ...V fcn C (x - y) J[ [5(x, - x) - 5(x, - y)] . 

(B9) 

Each leg of a vertex is labeled by a spatial variable Xj and an index fcj. The first two vertices 
computed from Eq. ( |B9|) are 



(x) = J(x) | J(y') V k C (x - y') dy', (BIO) 
F fc f(x,y) =-J(x)J(y)V fc V i C (x-y) + 5(x-y)J(x) ^ J(y')V fc V z C (x - y') dy' (Bll) 

(in the last expressions all integrations are made explicit, and we made use of the symmetry 
of the correlation function Co). Summation over indices and integration over spatial variables 
are performed when connecting lines, as usual. 



Since we are expanding in C g , the order of the diagram is equal to the number of the 
propagators. One can also obtain all diagrams of a given order by selecting the term of the 
given order in C g out of the expansion of Z g [H] in Eq. ( |B7p and acting with the functional 
derivatives on that term. Since each n-th power of C g is accompanied by the (2n)-th power 
of H, the only surviving derivatives will be of order not greater than 2n, which limits the 
possible vertices in the diagrams to a finite set. 

Since each vertex is quadratic in J, a connected diagram with n vertices V^ kl \ V^ kn ^ 
is of order (C g ) 1+ - +fcn anc l contributes to the cumulant G% n . The lowest-order diagrams are 
those with one propagator (Fig. lb); they yield leading-order expressions for C% = C and 
C4. Computing the diagrams in Fig. lb only requires the vertices and V^ 2 ' and gives 
Eqs. (|14|)-(|15|) after accounting for the factorial coefficients in Eq. ( \KQ ) and combinatorial 
factors. 



APPENDIX C: TWO-POINT CUMULANTS OF THE LENSED SKY 

The two-point cumulants Xmn can be straightforwardly obtained from the generating 
functional lnZ[J] of Eq. fljjD ; as described in Appendix [A], one needs to substitute the 
indicator function 

J pq (x)= P 8(x) + q6(x-9) (CI) 

into the generating functional and collect terms of the form p m q n . However, it is somewhat 
simpler to use Eq. (|C1|) directly with Eq. ([TT]) to derive the diagram rules specifically for 
the generating function Z x (p, q) of two-point cumulants, 

r 2 



z x (p> q) = { ex P 



Y(p 2 + q 2 )-pqCo[0 + g(9)-g(0)} 



(C2) 



Since only the difference Ag = g (9) — g (0) enters Eq. (|C2[ ) , we can simplify it by considering 
Ag as a random (vector-valued) variable whose distribution is specified by the generating 
function 
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Z A9 (h) =<expHh.(g(0)-g(O))]> 

— Z g [H — (5 (x — 0) — 5 (x)) h] = exp [~h 2 (C a (0) - C g (6))] . (C3) 

By a procedure similar to that we used above, we arrive to the following diagram rules. 
For each n — 1, 2, ... there is a vertex with n legs labeled by indices k± to k n , and the 
coefficient corresponding to the vertex is 



V n = pqV kl ...V kn C o (0), 



(C4) 



where the gradients are taken in the 2-dimensional Euclidean space, for instance 



V k C (x) = C (x) 



Xk 



w w r I \ n" I \ x k x l . f \ x $kl x k x l 

V k ViC (x) = C (x) — T + C (x) - 



x 



(C5) 
(C6) 



(Since we only consider small patches of the sky, we use Cartesian 2-dimensional vectors x 
interchangeably with the angular coordinate to denote points.) The "propagator" is 



2 [C g (0) - C g (9)] S kl = 2a] (0) 5, 



ki 



(C7) 



where the 5-tensor effectively sums over indices labeling the vertex legs. As usual in pertur- 
bation theory, one should count all diagrams differing by the way the vertices are connected 
and divide by the resulting combinatorial factor. 

The cumulant Xmn is obtained by adding all diagrams which give p m q n (with appropriate 
combinatorial factors) and multiplying by (m\n\). Since p and q enter the diagrams only 
in the combination pq, it follows that only cumulants Xmn with m = n are nonzero. The 
cumulant Xn is given by the sum of diagrams with one vertex; the lowest-order ones are the 
first diagram in Fig. 1(b) and the diagram in Fig. 2(a). These two diagrams give 



Xn(0) =C o (0) + a 
1 



a' 



2 1 1 

riii" i r' n r*n _i_ n' 



(C8) 



20 



(Here and below we imply that a 2 and all derivatives of the correlation function Cq (x) are 
evaluated at x = 6.) Diagrams with two vertices in Fig. 1(c) and Fig. 2(b) contribute to the 
cumulant X22 (0), 



Kf - ( f 



+ O {a e g ) . (C9) 



The diagram in Fig. 2(c) gives the leading-order expression for the cumulant X33, 



X33 (6) = 72a 4 g (C' ) 2 C' ( ; + O{a e g ). 



(CIO) 



Finally, estimate the variance of the cumulant X22 (#) for Gaussian skies. Because the 
temperature values are correlated, the set of all pairs of points separated by the angular 
distance 9 is not a statistically independent sample. The cumulant estimator X22 f° r a 
random function /(x) is a random variable expressed through /(x) as (cf. Eq. ([I] 



X22 



1 

A 



/ 2 (x)/ 2 (x')dxdx'-2 



/(x)/(x')dxdx' 



x-x'|= 



/ 2 (x)rfx 



(Cll) 



(here A is the area of the integration region). The expectation value of (X22) 2 involves at 
least a double integration over the observed region, 

<(X 22 ) 2 > = i / rfxrfx' / dydy'f 2 (x)f 2 (x')f(y)f(y') + ... (C12) 
A J|x-x'|=e "My-y'|=0 

The angle is typically much smaller than the separation between x and y in the integrals 

of Eq. ( |(J12j ), and we can replace x' by x, y' by y and so on in that equation: 



;(X2 2 ) 2 > / 4 (x)/ 4 ( y )dxrfy-A J / 4 (x)/ 2 (y)/ 2 (z)dxrfydz 

+ J f 2 (w)f 2 Wf 2 (y)f 2 ( Z )dwd X dydz, 



(C13) 
(C14) 



where the integrations in w, x, y, z are over all sky. This approximation is equivalent to 
replacing X22 (#) by the one-point cumulant 



X4 = </ 4 >-3</7 
21 



(C15) 



which is adequate for estimating the approximately ^-independent variance of X22 (9) at 
small 9. If the distribution of / is Gaussian with zero mean and correlation function Co (x), 
then Eq. ( |C14 ) gives after some algebra 



<(X22) 2 ) 



2 \ 432 
A 3 



36 



J C (x) C (y) Co (z - x) C (z - y) rfxdydz 

24 r 

+ ^ y C$ (x) dx. (C16) 



C 2 (x) cix 



The first two integrals can be expressed through the angular power spectrum, and we obtain 

2 



{(X22) 2 ) 



432 



36 



5>/ + l)C? 



/»7T 

+ 12 / C%(0)&n9dB. (C17) 



This is the final formula we use for the variance of the cumulant estimator X22 {9) at small 
angles 9. 

APPENDIX D: AN ALGORITHM FOR COMPUTING iV-POINT CUMULANTS 



Cumulants of a distribution of n variables labeled by n indices fci, fc„ 

and can be defined by 



X.k\...k n — I * 



d_ 

dpi 



. d 

dp n 



k/1 



ln/(pi, ...,p r , 



(Dl) 



where / is the generating function for the joint distribution of Xj. We call k\ + ... + k n the 
order of the cumulant Xk 1 ...k n - The moments of the distribution are expressed through the 
same generating function as 



i-kn — ( X l 1 --- X n ra ) — ( * 



<9 



fcl 



<9 



f(pi,...,Pn) ■ 



(D2) 



a%=0 



Although Eqs. ( |U1| ) ( [02]) allow to express all cumulants through the moments, a direct 
calculation for all combinations of indices is cumbersome to implement. We present a com- 
pact recursive algorithm for computing all x' s whereby the cumulants corresponding to all 
combinations of k's are evaluated using the previously computed cumulants of lower orders. 
For the case n = 1, this algorithm is given in Ref. ]16| and is based on the identity 
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»s = Yl (l _ i)xkHs-k- (D3) 

fc=0 ^ ' 

Here we imply that xo = and //q = 1, and that binomial coefficients with negative entries 
vanish. From Eq. (|D3|) we can recursively express Xs through the lower-order cumulants, 

5-1 

fc=0 

In the general case n > 1, we need an analog of the identity of Eq. ( |D3j ). The resulting 
formula turns out to be 

Ev^ ki + ... + k n f sA ( s n \ /T\t:\ 
■■■ 7. : : I, )■■■[ , Xfci...fcn/ i si-fci,..., s „-fc„, {Mo) 

ki=0 k n =0 x ' x ' 

from which one obtains the cumulant Xs\...s n as a function of moments \x and cumulants 



Y^(l-i) Xk,i '- k - (D4) 



X of lower orders. Using Eq. ( p5|) , one can organize the computation of all cumulants 
order by order in an iterative manner. Given iV sample realizations of the set of random 
variables {xi}, one first calculates the moments fi Sl ... Sn up to the required order; this takes 
a number of operations linear in N. The first-order cumulants Xioo...o> X010...0 etc. are equal 
to the similarly labeled moments. Thus having obtained the base of induction, one iterates 
through increasing order of the cumulants index by index, using Eq. ( P5|) to express the 
next cumulant through moments and previously computed cumulants of lower orders. 

Now we outline the derivation of Eq. ( [D 5| ) . We combine Eqs. ( |D1|) and (|D2| ) into an 
identity 



cxp 



00 ki h 

2^Xk x ...k n k , ± , 

M=0 



Jti ^k„ 



ki=0 1 n [k] [ 1 



Regarding Hk x ...kn = A*[fc] as functions of x's, we take the partial derivative of Eq. ( |D6| ) with 
respect to X[k] f° r a fixed multi- index [k]. By equating the coefficients at p{ 1 ...p^ > = [p s ] we 
obtain the relation 

d [si! 



dxik] m W[s-k]V 

where the multi-index [s — k] means {(si — k) , (s„ — A; n )} in our condensed notation. We 
also introduce for convenience the generalized binomial coefficient 
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t*n_ w /v, (D8) 



JA;]7 [jfe]![s-jfe]! VW \*« 
The form of Eq. ( P3|) leads us to conjecture that the generalized identity should look 

like 

fc 1= fc n =o ^ ' 

with some combinatorial coefficients {^} which should be similar to those of Eq. ( |D8| ). The 
requirement that these coefficients must satisfy Eq. ( [D 7|) leads after some algebra to 

£k|L-M (1)1(1) 



which yields Eq. (|D5|) after substitution into Eq. ( P9|) . 



24 



REFERENCES 

[1] G. Jungman, M. Kamionkowski, A. Kosowsky, D. N. Spergel, Phys. Rev. D 54, 1332 
(1996). 

[2] C. L. Bennett, Bull. Am. Astron. Soc. 189, 8805 (1996). See also the MAP page "http: 
/ /map. gsfc.nasa.gov/" . 

[3] The PLANCK project, "http: //astro. estec.esa.nl /SA-general /Projects /Planck/". 

[4] U. Seljak, Ap. J. 463, 1 (1996). 

[5] E. V. Linder, Astron. Astrophys. 323, 305 (1997); M. Suginohara, T. Suginohara, and 
D. N. Spergel, Ap. J. 495, 511 (1998). 



[6] R. Stompor and G. Efstathiou, preprint |astro-ph/ 9805294] (1998 



[7] F. Bernardeau, Astron. Astrophys. 324, 1 (1997). 
F. Bernardeau, preprint |astro-ph/ 98022^3] (1998). 



[9] P. G. Ferreira, J. Magueijo, and J. Silk, Phys. Rev. D 56, 4592 (1997). 
[10] S. Winitzki and A. Kosowsky, New Astronomy 3, 75 (1997). 

[11] J. Schmalzing and T. Buchert, Ap. J. 482, LI (1997); J. Schmalzing and K. M. Gorski, 
preprint |astro-ph/97lOT85| (MNRAS, in press, 1998). 



[12] U. Seljak and M. Zaldarriaga, Ap. J. 469, 437 (1996); M. Zaldarriaga and U. Seljak, 
preprint |astro-ph /9803 1 50| (1998). The CMBFAST program is available at "http: //arc- 



turus.mit.edu:80 /~matiasz /CMBFAST /cmbfast.html" . 
[13] U. Seljak, Ap. J. 436, 509 (1994). 

[14] X. Luo and D. Schramm, Phys. Rev. Lett. 71, 1124 (1993); A. Gangui, F. Lucchin, 
S. Mattarese, and S. Mollerach, Ap. J. 430, 447 (1994). 



25 



[15] P. Coles, MNRAS 234, 509 (1989); J. R. Gott et al, Ap. J. 352, 1 (1990); G. Smoot 
et al, Ap. J. 437, 1 (1994). 

[16] M. G. Kendall and A. Stuart, The Advanced Theory of Statistics, 4th ed., Charles Griffin 
(1977). 



26 



FIGURES 




(a) (b) 

FIG. 1. Diagram representation of the generating function for the lensed sky. (a) A vertex 
labeled n has n legs, corresponding to ra-th derivative of Co- The propagator corresponds to the 
lensing correlation function C g . (b) Lowest-order diagrams with one propagator lead to corrections 
in the power spectrum and a non-Gaussian 4-point correlation function. 




(a) (b) (c) 

FIG. 2. Diagrams leading to second-order effects: (a) correction to the power spectrum, (b) 
second-order terms in the cumulant X22, (c) leading-order term in the cumulant X33- 



n 



1 v- 

+2 

1 \ x 



(a) 



(b) 



FIG. 3. Diagrams for non-Gaussian initial conditions, (a) In addition to previous types of 
diagrams, there are "propagators" with more than two legs which correspond to non-Gaussian 
lensing sources, (b) An example of a new diagram contributing to the 4-point cumulant. 
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0.1 1 10 

Angular scale (degrees) 

FIG. 4. The dimensionless cumulants \22 {&) and X33 {&) as functions of the angular scale 6. 
The cosmological parameters are: for solid lines Qf, = 0.05, Qa = 0.3, h = 0.5, n s = 1.2 (open); 
for dashed lines Vt^ = 0.03, Oa = 0.97, h = 0.7, n s = 1.0 (flat). The horizontal lines represent the 
standard deviation of \22- 
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